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SUMMARY 

The  results  of  an  experimental  study  aimed  at  characterizing  the  behavior  of  concrete  for  high  confining 
pressures  (up  to  500  MPa)  are  reported.  The  main  characteristics  of  the  response  under  deviatoric  conditions 
are  quasi-linearity  in  the  elastic  regime,  stress-path  dependency,  and  gradual  change  from  compressibility 
to  dilatancy  under  increasing  deviatoric  stress.  By  performing  the  cyclic  triaxial  compression  tests  with 
several  load-creep-unload  and  reload  cycles,  the  time  influence  on  the  overall  behavior  was  detected. 
Further,  whether  the  main  features  of  the  observed  behavior  can  be  described  within  the  framework 
of  elastic/viscoplasticity  theory  was  investigated.  For  this  purpose,  Cristescu’s  ( Rock  Rheology.  Kluwer 
Academic  Publishers:  The  Netherlands,  1989)  approach  was  used.  It  was  shown  that  the  proposed 
elastic/viscoplastic  model  captures  the  main  features  of  concrete  behavior  at  high  pressures.  Copyright 
©  2008  John  Wiley  &  Sons,  Ltd. 
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1.  INTRODUCTION 

A  rather  significant  literature  exists  on  the  quasi-static  mechanical  response  of  cementitious  mate¬ 
rials  and  geomaterials.  However,  most  of  the  experimental  results  are  obtained  using  triaxial  cells 
and  concern  the  response  under  moderate  confining  pressures,  i.e.  up  to  approximately  100  MPa 
(see,  for  example,  [1-4],  etc.).  True  triaxial  compression  test  results  reported  are  for  the  same  range 
of  pressures  (see,  for  example,  [4, 5]).  Bazant  et  al.  [6]  have  developed  an  oedometric-type  test  that 
allows  characterization  of  the  response  at  high  pressures.  A  cylindrical  specimen  (height  44.4  mm, 
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diameter  18.8  mm)  is  placed  in  a  very  rigid  hollow  plate,  which  prevents  lateral  deformation,  and 
is  subjected  to  compression  parallel  to  its  vertical  axis.  The  major  limitations  are  related  to  the  size 
of  the  specimens  and  accurate  evaluation  of  frictional  effects.  Bourlion  et  al.  [7]  have  extended 
this  testing  method  to  allow  testing  of  concrete  with  aggregate  sizes  up  to  16  mm  at  levels  of  mean 
pressure  of  about  600  MPa.  These  authors  have  also  reported  true  hydrostatic  test  results  on  small 
mortar  size  specimens. 

There  is  a  large  body  of  literature  on  modeling  the  behavior  of  cementitious  materials.  For 
moderate  levels  of  pressure,  the  non-linearity  of  the  response  is  mainly  due  to  damage  by  micro- 
and  meso-crack’s  growth.  Major  developments  have  occurred  over  the  past  decades  in  modeling 
the  behavior  at  continuum  level.  Continuum  damage  models  of  various  degrees  of  complexity 
have  been  proposed  (for  comprehensive  reviews  of  the  main  contributions,  see  [8,  9]).  For  higher 
levels  of  pressure  (mean  stress  or  confinement),  the  behavior  is  more  ductile  and  plasticity 
models  (e.g.  [10,11])  as  well  as  coupled  damage-plasticity  models  ([12-16],  etc.)  are  more 
appropriate. 

To  model  time  effects  on  the  non-linear  behavior  of  cementitious  materials,  a  very  large 
number  of  models  developed  for  inviscid  behavior  have  been  extended  in  the  rate-dependent 
range  using  the  overstress  concept  (see,  [17]).  Thus,  the  expressions  of  the  constitutive  func¬ 
tions  (i.e.  yield  function,  plastic  potential)  used  in  these  rate- sensitive  models  are  the  same  as 
the  ones  used  to  describe  inviscid  behavior.  In  contrast,  the  approach  introduced  by  Cristescu 
(see,  for  instance,  [4])  is  to  determine  the  specific  expressions  of  the  yield  function  and 
viscoplastic  potential  for  a  given  material  directly  from  creep  data.  This  allows  describing  both 
compressibility  and  dilatancy  by  a  unique  continuous  yield  function  and  viscoplastic  potential, 
respectively. 

In  this  paper,  we  present  the  results  of  an  experimental  investigation  into  the  effects  of  high 
confinement  on  the  deformation  and  strength  of  a  concrete  material.  By  performing  cyclic  triaxial 
compression  (CTC)  tests  with  several  load-creep-unload  and  reload  cycles,  the  time  influence 
on  the  overall  behavior  was  detected.  Based  on  these  data,  an  elastic/ viscoplastic  model  that 
captures  compressibility,  dilatancy,  and  strain  rate  effects  has  been  developed.  The  structure  of  this 
model  and  a  procedure  for  determining  the  expressions  of  the  constitutive  functions  and  material 
parameters  based  on  a  minimal  set  of  data  are  given  in  Section  3.  We  conclude  by  presenting 
comparison  between  model  simulations  and  data. 


2.  EXPERIMENTAL  STUDY 

The  aim  of  this  experimental  study  is  to  extend  the  existing  database  on  the  deformation  and 
failure  of  concrete  by  investigating  the  behavior  over  a  range  of  pressures  characteristic  of  impact 
events.  Of  particular  interest  was  the  influence  of  very  high  confining  pressure  (up  to  500  MPa)  on 
the  material’s  deformation  and  strength.  Thus,  the  experimental  program  included  (a)  quasi-static 
hydrostatic  compression  tests  from  which  the  increase  in  the  compressibility  of  the  material  and 
an  estimate  of  the  compaction  properties  over  the  range  pressures  encountered  in  impact  events 
can  be  inferred,  (b)  quasi-static  CTC  tests  under  lateral  confining  pressure,  <73,  ranging  from  50 
to  500  MPa,  based  on  which  the  effect  of  loading  history  on  the  response  and  the  importance  of 
short-term  time  effects  on  the  behavior  can  be  determined.  A  total  of  21  tests  were  conducted 
at  the  U.S.  Army  Corps  of  Engineer  Engineering  Research  and  Development  Center,  Waterways 
Experiment  Station  (WES)  in  Vicksburg,  Mississippi. 
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2.7.  Material  studied 

The  mixture  proportions  and  resulting  basic  properties  of  the  concrete  studied  are  given  in  Table  I. 
The  concrete  is  a  mixture  of  Portland  cement,  chert  aggregates  (fine  and  coarse),  fly  ash,  water,  and 
water-reducing  admixture.  The  uniaxial  compressive  and  tension  strengths  at  28  days  were  45.6 
and  3. 13  MPa,  respectively.  The  test  results  indicate  that  under  uniaxial  compression,  the  material 
exhibits  typical  brittle  behavior. 

2.2.  Specimen  preparation  and  quasi- static  test  instrumentation 

Prior  to  each  test,  the  prepared  specimen  was  weighed,  carefully  measured,  and  ultrasonic  measure¬ 
ments  (P-  and  S-wave  velocities)  were  taken.  For  all  tests,  the  specimens  had  a  nominal  height  of 
110  mm  and  a  diameter  of  50  mm.  Because  of  the  inherent  heterogeneity  of  the  material  and  the 
high-pressure  range  in  which  it  was  tested,  instrumentation  with  strain  gages  was  found  unreli¬ 
able.  Instead,  to  measure  the  displacement  between  the  top  and  the  base  caps  (i.e.  determine  the 
axial  deformation  of  the  specimen),  two  linear  variable  differential  transformers  (LVDTs)  mounted 
vertically  on  the  instrumentation  stand  and  positioned  180°  apart  were  used.  The  radial  deflection 
was  measured  using  two  LVDTs,  which  were  connected  to  two  small  steel  footings  mounted 
diametrically  (180°  apart)  at  the  mid-height  of  the  specimen  (see  Figure  1).  A  gypsum  paste  was 
used  to  patch  all  voids.  Subsequent  to  patching  of  the  voids,  the  specimen  was  encased  in  a  latex 
membrane,  and  this  membrane  was  then  surrounded  by  a  layer  of  putty-like  material.  The  putty 
served  the  function  of  self- sealing  any  potential  punctures  of  the  inner  membrane.  A  second  latex 
membrane  was  then  applied  over  the  putty  and  the  membrane  was  coated  with  a  sealant  to  prevent 
its  degradation  by  the  petroleum-based  hydraulic  fluid. 

The  axial  load  was  provided  by  an  8.9  MN  loading  machine  (Secodyne)  coupled  with  600- 
MPa-capacity  pressure  cell  (see  the  schematic  diagram  in  Figure  2).  The  loading  rate/profile  was 
computer  programmable,  allowing  for  precise  control  of  the  experiments  under  load,  displacement, 
or  strain  control. 

2.3.  Experimental  results 

A  triaxial  compression  test  is  generally  conducted  in  two  phases:  a  hydrostatic  compression 
phase  and  a  deviatoric  phase.  During  the  hydrostatic  phase,  the  specimen  is  subjected  to  isotropic 
compression,  while  measurements  of  the  axial  and  radial  deformations  are  made.  The  deviatoric 


Table  I.  Ingredients  and  mixture  proportions  for  WES 5 000  concrete. 


Item 

Mixture  proportions, 
saturated  surface-dry 

Type  I  Portland  cement 

264.0kg/m3 

Flyash 

55.8  kg/m3 

9.5  mm  local  unprocessed  chert  coarse  aggregate 

1037.6kg/m3 

Local  unprocessed  chert  fine  aggregate 

840.7  kg/m3 

Water 

145.9kg/m3 

Water-reducing  admixture  ‘300N’ 

0.651/m3 

High-range  water  reducing  admixture  ‘Rheobuild  716’ 

1.61/m3 
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Figure  1.  Spring-arm  lateral  LVDT  mounted  on  test  specimen. 


phase  of  the  test  is  conducted  after  the  desired  confining  pressure  has  been  reached  during  the 
hydrostatic  phase.  While  holding  the  desired  lateral  confining  pressure  constant,  the  axial  load  is 
increased  until  the  specimen  fails.  As  an  example,  the  pressure  vs  volumetric  data  obtained  from 
the  hydrostatic  phase  of  a  CTC  test  under  200  MPa  confining  pressure  are  presented  in  Figure  3.  In 
order  to  evaluate  the  bulk  modulus,  an  unloading-reloading  cycle  was  performed  at  100  MPa.  Note 
the  significant  hysteresis  displayed,  which  clearly  makes  the  determination  of  the  bulk  modulus 
highly  imprecise. 

To  determine  the  elastic  properties  with  accuracy,  it  is  necessary  to  ensure  a  good  separation 
between  viscous  effects  and  unloading.  One  way  of  accomplishing  this  is  to  hold  the  stress  constant 
before  each  unloading/reloading  cycle  such  as  to  allow  the  material  to  creep.  The  duration  of  the 
creep  stage  (time  interval  in  which  the  load  is  kept  constant)  is  usually  determined  by  monitoring 
the  rate  of  change  of  the  strains.  When  the  strain  rate  approaches  zero,  i.e.  the  material  reaches 
by  creep  a  quasi-stable  state,  unloading  is  conducted.  More  details  concerning  this  experimental 
procedure  and  its  application  to  a  variety  of  geologic  materials  can  be  found  in  Cristescu  and 
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Figure  2.  Section  through  triaxial  pressure  cell. 


Figure  3.  Mean  stress-volumetric  strain  data  from  a  hydrostatic  cyclic  test  on  concrete. 
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Hunsche  [4].  This  step-wise  loading  procedure  was  used  for  the  hydrostatic  phase  of  each  CTC  test 
conducted.  As  an  example,  the  results  of  a  hydrostatic  test  up  to  500  MPa  are  shown  in  Figure  4. 

The  close  agreement  between  the  axial  and  radial  gage  measurements  indicates  that  the  material 
is  isotropic.  The  test  consisted  of  seven  loading-creep-unloading-reloading  cycles.  For  each  cycle, 
prior  to  unloading,  the  pressure  was  held  constant  for  30  min.  Then,  partial  unloading  and  further 
reloading  were  performed.  In  Figure  5,  the  volumetric  strain  rate  vs  time  data  corresponding  to  the 


Strain  (mm/mm) 


Figure  4.  Mean  stress-volumetric  strain  data  from  a  hydrostatic  test  up  to  500  MPa  conducted  following 
a  step-wise  loading/creep/partial  unloading  procedure. 


2600  2600  3000  32TO  3400  3500  3500 

Time  (sec) 

Figure  5.  Volumetric  strain  rate  vs  time  from  the  second  cycle  (unload  at  100  MPa)  of  a 

hydrostatic  test  up  to  500  MPa. 
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second  cycle  (unloading  at  100  MPa)  are  presented.  It  can  be  seen  that  within  the  duration  of  the 
creep  stage,  the  rate  of  change  of  the  volumetric  strain  with  time  approached  zero.  In  Figure  6, 
the  pressure-strain  curve  corresponding  to  the  second  creep-loading-unloading-reloading  cycle 
is  shown.  Note  that  very  little  hysteresis  is  present,  and  the  slope  of  the  quasi-linear  part  of  the 
unloading  curve  coincides  with  the  slope  of  the  quasi-linear  part  of  the  reloading  curve.  The  bulk 
modulus,  K,  corresponding  to  several  pressure  levels  was  evaluated  from  the  slopes  of  the  quasi- 
linear  part  of  the  unloading  curves  of  pressure  vs  volumetric  strain.  The  bulk  modulus  was  found 
to  increase  with  the  applied  pressure  from  K  =  23.24GPa  (second  cycle)  to  K  =  28.76GPa  (last 
cycle).  The  reduction  in  porosity  at  the  end  of  the  test  (i.e.  at  0.5  GPa)  is  approximately  2.5%. 
Figure  7  shows  the  principal  stress-strain  data  from  the  deviatoric  phase  of  a  CTC  test  conducted 
at  (j3  =  375  MPa  confining  pressure  following  a  step-wise  loading  procedure.  Note  that  the  material 
exhibits  irreversible  time-dependent  behavior.  The  stress  vs  volumetric  strain  curve  also  shows  a 
highly  non-linear  volume  response,  with  the  onset  of  dilatancy  being  very  close  to  failure.  The 
same  trends  have  been  observed  in  all  the  tests  conducted.  Figure  8  shows  a  comparison  of  the 
stress-strain  data  from  the  deviatoric  phase  of  CTC  tests  at  confining  pressures  of  200,  300,  and 
450  MPa,  whereas  in  Figure  9  presents  the  strength  values  determined  in  the  same  tests.  It  is  clearly 
seen  that  higher  confinement  produces  higher  yield  stress,  increased  hardening,  and  increased 
strength. 

From  the  test  results,  we  can  conclude  that  for  high  pressures,  the  material  displays  irreversible 
time-dependent  properties.  On  the  other  hand,  the  analysis  of  the  variation  of  Young’s  modulus  with 
pressure  (see  Figure  10)  shows  that  for  specific  levels  of  radial  confinement,  E  is  decreasing 
with  the  mean  stress,  which  is  indicative  of  damage.  Modeling  of  damage  and  damage  couplings 
with  plasticity  and  viscous  effects  is  beyond  the  scope  of  this  paper.  We  will  model  only  the 
main  features  of  the  behavior,  i.e.  time  effects  on  the  plastic  behavior  within  the  framework  of 
viscoplasticity. 


Figure  6.  Mean  stress-volumetric  strain  data  from  the  second  cycle  (unload  at  100  MPa)  of  a 

hydrostatic  test  up  to  500  MPa. 
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Figure  7.  Stress-strain  data  from  the  deviatoric  phase  of  a  cyclic  triaxial  compression  test  at  a 

confining  pressure  of  375  MPa. 


Figure  8.  Stress-strain  data  from  the  deviatoric  phase  of  cyclic  triaxial  compression  tests  at  confining 

pressures  of  200,  300,  and  450  MPa. 


3.  ELASTIC-VISCOPLASTIC  MODEL  DEVELOPMENT 

The  most  widely  used  approach  to  model  coupling  between  plastic  and  viscous  effects  on  the 
behavior  of  geological  or  cementitious  materials  is  that  proposed  by  Perzyna  [17]  (see  also  [18]). 
The  basic  assumption  is  that  the  viscous  properties  of  materials  become  manifest  only  after  the 
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Figure  9.  Ultimate  strength  as  a  function  of  confining  pressure. 


Figure  10.  Comparison  between  theoretical  and  experimental  variations  of  Young’s 
modulus  E  with  the  mean  stress. 


passage  to  the  plastic  state.  The  expressions  of  the  constitutive  functions  (yield  function,  viscoplastic 
potential)  are  the  same  as  the  ones  used  to  describe  the  elastic/plastic,  time-independent  behavior. 
In  this  paper,  we  will  use  the  elastic/viscoplastic  modeling  framework  proposed  by  Cristescu 
(see,  for  instance,  [4]),  i.e.  the  mathematical  expressions  of  the  constitutive  functions  are  a  priori 
unknown  and  will  be  determined  from  the  experimental  data  (creep  test  results).  We  begin  by  briefly 
reviewing  the  general  elastic/viscoplastic  framework  and  then  present  the  simplifying  assumptions 
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made  in  order  to  obtain  a  reasonably  simple  mathematical  model  for  the  given  material.  It  is 
assumed  that  the  material  is  homogeneous  and  isotropic.  The  reference  configuration  is  the  actual 
configuration.  The  displacements  and  material  rotations  are  assumed  to  be  small.  Thus,  the  strain 
rate  £  can  be  decomposed  additively  into  an  elastic  part,  £E,  and  an  irreversible  part,  a1,  i.e. 

£  =  £E  +  £J  (1) 


The  elastic  component  of  the  strain  rate  is  given  by 


£E 


_G_  (\ _ l_\ 

~  2G  ~\3/<r  2  g)P 


(2) 


where  the  dot  stands  for  the  derivative  with  respect  to  time,  G  and  K  are  the  shear  and  bulk  moduli, 
respectively,  I  is  the  second-order  identity  tensor,  and  p  is  the  mean  stress  (i.e.  p  =  tr(a)/3).  No 
assumption  concerning  the  existence  of  a  viscoplastic  potential  is  made.  Thus,  the  irreversible 
strain  rate  is  considered  to  be  given  by  an  overstress-type  law  of  the  form: 


&r(l  — 


W(t) 

H(<j) 


O) 


(3) 


where  t  is  the  actual  time  and  W(t)  is  the  irreversible  stress  work  per  unit  volume.  The  only 
restriction  imposed  on  the  tensor- valued  N(g),  which  defines  the  orientation  of  the  viscoplastic 
strain  rate,  is  to  be  isotropic  (see  [19]).  In  (3),  the  symbol  (),  known  as  the  Macaulay  bracket, 
denotes  (x)  =  (x  +  \x\)/2  and  kj  is  a  viscosity  coefficient.  Thus,  there  is  viscoplastic  flow  if  and 
only  if  H(a)>W(t).  Only  transient  creep  is  modeled.  From  (3)  it  follows  that  the  deformation 
due  to  transient  creep  stops  after  a  certain  finite  time  interval  and  stabilization  takes  place,  the 
equation  for  the  stabilization  boundary  for  creep  being 


H(g)  =  W(t) 


(4) 


Note  that  this  choice  of  hardening  variable  is  appropriate  for  the  description  of  the  behavior  of 
concrete  because  it  captures  both  the  shear  and  the  volumetric  inelastic  effects.  Indeed,  W(t)  can 
be  decomposed  into  two  terms: 

W(t)=  f  a(t):&(t)dt  =  WH(t)  +  WD(t)  [  p(t)^(t)dt+ f  a'(t):i\t)dt  (5) 

Jo  Jo  Jo 


where  Wu(t)  is  the  energy  related  to  irreversible  volume  change  that  is  stored  during  compaction  or 
released  during  dilatancy  (volumetric  expansion)  of  the  material,  and  is  always  positive  and 

represents  the  energy  input  needed  for  change  in  shape.  In  (5),  £*  (0  is  the  irreversible  volumetric 
strain  rate  and  a\t)  is  the  stress  deviator. 

The  first  step  in  the  development  of  the  elastic-viscoplastic  model  for  the  material  studied  is  the 
description  of  the  elastic  behavior.  For  each  CTC  test  conducted,  Young’s  modulus,  E ,  and  the  bulk 
modulus,  Ky  were  evaluated  from  the  slopes  of  the  quasi-linear  parts  of  the  unloading-reloading 
cycles  performed  at  the  end  of  each  of  the  creep  steps.  We  also  estimated  the  corresponding  Poisson 
ratio,  v.  As  it  is  the  case  with  most  rocks,  soils,  and  cementitious  materials,  the  elastic  moduli  are 
stress  dependent.  Thus,  the  stress-strain  relationship  (2)  is  hypoelastic  in  the  sense  of  Truesdell 
[20].  The  necessary  and  sufficient  conditions  for  an  isotropic  material  with  stress-dependent  elastic 
moduli  to  have  a  conservative  response  in  the  elastic  regime  were  established  by  Loret  [21].  It  was 
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shown  that  the  choice  of  the  laws  of  variation  of  the  elastic  moduli  with  the  stress  state  cannot  be 
arbitrary.  The  laws  of  variation  should  be  expressible  in  terms  of  the  stress  invariants:  the  mean 

stress  p ,  and  the  second  invariant  of  the  stress  deviator  q  =  -J jtr(a')2 .  An  analysis  of  the  variation 
in  Young’s  modulus  with  pressure  indicates  a  global  increase  in  E  with  increasing  p.  However,  for 
certain  levels  of  radial  confinement,  a  decrease  in  E ,  which  is  indicative  of  damage,  is  observed. 
Modeling  of  damage  and  damage  couplings  with  plasticity  and  viscous  effects  is  beyond  the  scope 
of  this  paper.  We  model  only  the  increase  in  E  with  pressure  by  considering  the  following  law  of 
variation: 


E  (p)  =  E°°  —  b  exp  i-p/c)  (6) 

In  (6),  E°°  is  the  asymptotic  limit  towards  which  E  tends  as  the  pressure  becomes  very  large, 
(p  — ►  oo ),  i.e.  the  value  of  the  Young  modulus  corresponding  to  a  fully  compacted  state.  Good 
overall  agreement  with  all  the  available  data  is  obtained  for  £00  =  61GPa,  b  =  E°°  —  E\p=o  = 
18.92GPa,  and  c  =  b(dE/dp\p=o)  =  0.33 MPa  (see  Figure  10).  We  consider  v  =  0.175,  a  value  that 
corresponds  to  the  mean  of  the  experimental  values  from  all  tests  (see  Figure  11),  and  calculate 
K  from  the  elasticity  relationship: 


K  = 


E 

3(1 -2v) 


(7) 


Several  factors  entered  into  assuming  a  constant  value  of  the  Poisson  ratio.  First,  for  the  sake 
of  simplicity,  it  was  seen  advisable  to  limit  the  number  of  coefficients  involved  in  the  model. 
Also,  due  to  the  overall  higher  level  of  confidence  in  the  axial  measurements,  it  was  found  best  to 
approximate  the  experimental  variation  of  E  with  the  mean  stress  rather  than  to  approximate  the 
experimental  evolution  of  the  bulk  modulus.  Furthermore,  a  value  of  0.175  for  the  Poisson  ratio  is 
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Figure  11.  Experimental  variation  of  Poisson’s  ratio  with  the  mean  stress. 
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well  within  the  range  of  values  typically  observed  for  concrete  (see  [22]).  Note  that  Equation  (7) 
with  E  given  by  (6)  approximates  well  the  experimental  variation  of  the  bulk  modulus  K  with 
the  mean  stress  (see  Figure  12).  The  material  being  isotropic,  the  yield  function,  77,  depends  on  a 
through  its  invariants.  The  experimental  data  available  were  obtained  from  CTC  tests  in  which  only 
two  principal  stresses  can  be  varied  and  measured  independently,  i.e.  o\  and  02  =  o?>-  The  effect 
of  the  third  invariant  of  stress  cannot  be  inferred  without  data  for  which  all  three  principal  stresses 

can  be  varied  independently.  That  is  why  we  assume  that  H(a)  =  H(p,q),  where  q  =  -yj^ tr(a7)2 
is  the  second  invariant  of  the  stress  deviator  g'.  To  determine  the  specific  expression  of  H(p,q) 
based  on  CTC  results,  we  adopt  the  procedure  proposed  by  Cristescu  [3]  (see  also  [4]).  Thus,  we 
assume  that  H(p,q)  can  be  expressed  as 

H(p,q)  =  HH(p)  +  HD(p,q)  (8) 

such  that 

HD(p,  0)=0  (9) 

where  Hu(p)  =  Wu  and  H&(p,  q)  =  Wb,  Wn  and  Wd  being  the  volumetric  and  deviatoric  parts  of 
the  irreversible  stress  work  at  creep  stabilization.  Thus,  in  order  to  determine  Hu(p ),  we  calculated 
the  irreversible  stress  work  corresponding  to  creep  stabilization  according  to  the  following  equation: 

wh«h)=  r  pmliodt  do) 

Jo 

using  data  from  the  hydrostatic  phase  of  CTC  tests  at  150,  300,  450,  and  500  MPa  confinement, 
respectively.  In  (10),  s\  is  the  irreversible  volumetric  strain  and  tn  corresponds  to  the  end  of  the 
hydrostatic  phase  of  a  CTC  test.  Next,  we  plotted  the  obtained  values  of  Wn(t)  as  a  function 
of  p.  Note  that  Wu  increases  monotonically  with  the  mean  pressure.  However,  for  very  high  levels 


Figure  12.  Comparison  between  theoretical  and  experimental  variations  of  the  bulk 
modulus  K  with  the  mean  stress. 
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of  pressures,  the  Wu  vs  p  curve  should  reach  a  plateau,  which  corresponds  to  a  state  where  all 
the  micro-cracks  and  pores  are  closed  (i.e.  material  is  fully  compacted).  This  plateau  obviously 
lies  beyond  the  hydrostatic  pressure  limit  of  the  equipment  available.  An  approximate  value  may 
be  inferred  from  post-test  observations  of  the  density  distribution  of  concrete  targets  impacted  at 
conventional  striking  velocities.  Indeed,  X-ray  measurements  show  that  at  about  2-3  projectile 
radii  around  the  penetration  tunnel,  the  target  material  is  fully  compacted.  Assuming  that  the 
pressure  corresponding  to  fully  compacted  state  is  of  the  same  order  of  magnitude  as  the  penetrator 
strength,  the  plateau  is  at  about  1  GPa  (see  [23]).  The  experimental  Wu  data  at  stabilization  were 
approximated  with  the  function: 


Hu(p)  = 


(11) 


with  ah  =  12.6MPa,  bh  =364.28 MPa,  Ch  =  —  365MPa,  *4  =30.7  MPa,  and  eh  =0.003 MPa. 
Figure  13  shows  the  values  of  the  irreversible  volumetric  stress  work  obtained  from  experimental 
data  in  the  hydrostatic  tests  and  the  approximating  function. 

Similarly,  Hv(p,q)  is  determined  from  the  data  obtained  in  the  deviatoric  phase  of  the  CTC 
tests.  The  irreversible  stress  work  is  computed  using 


Wd(T)=  f  <73(04(0*+  f  t^tH  (12) 

Jtu  Jtn 

where  tu  represents  the  beginning  of  the  deviatoric  part  of  the  test  and  T  corresponds  to  the  end 
of  the  creep  stage.  Figure  14  shows  the  experimental  values  of  Wd  obtained  from  several  CTC 
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Figure  13.  Stabilization  boundary  as  a  function  of  the  mean  stress  (solid  line);  experimental  values  of  the 
irreversible  stress  work  at  creep  stabilization  for  various  confining  pressures  (symbols). 
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Figure  14.  Experimental  values  of  the  irreversible  stress  work  at  stabilization  as  calculated  based  on  the 
results  in  the  deviatoric  phase  of  several  CTC  tests. 


tests.  At  any  given  confining  pressure,  the  data  trend  is  the  same,  i.e.  initially  linear,  eventually 
curving  upward,  and  asymptoting  towards  infinity  as  the  specimen  approaches  failure. 

A  function  that  well  approximates  the  data  is 

HD(p,q)=d0q+diq2+d2exp(d3q2)-d2  (13) 

The  coefficients  d,  depend  on  the  confining  pressure  (T3  and  are  expressed  as 


do  (03)  =  d0a  ■  (1  -exp(<%<73)) 
d\(o3)  =d\a-(\- exp(d\ho3)) 
d2((J3)  =  d2a  +d2b  ■  (1  -exp(d2ca3)) 


d3(<73)  = 


\+d3-d3,-a3 


(14) 

(15) 

(16) 

(17) 


with  doa  =0.0141,  dob  =  —0.0075,  d\a  =4.08 x  10_5MPa_1,  dh  =  -0.005 MPa-1,  d2a  =  1.6x 
10_4MPa,  d2h  =0.0564 MPa,  d2c  =  —0.0082 MPa- 1 ,  d3a  =0.00 135 MPa-1,  and  d3b  =  175.26.  As 
an  example,  a  comparison  between  the  approximating  function  and  the  experimental  results 
corresponding  to  200  MPa  confining  pressure  is  shown  in  Figure  15.  The  results  of  this  test  were 
not  used  in  the  determination  of  (13).  A  very  good  agreement  is  observed. 

According  to  the  constitutive  equation  (3),  the  tensorial  function  N{g)  governs  the  orientation 
of  the  irreversible  strain  rate,  e1.  For  an  isotropic  material,  N(g)  must  satisfy  the  invariance 
requirement 

N(QaQT)  =  QN(a)QT  (18) 
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Figure  15.  Comparison  between  theoretical  and  experimental  variations  of  the  irreversible  stress  work  as 
a  function  of  the  deviatoric  stress  at  a  confining  pressure  of  200  MPa. 


for  any  orthogonal  transformation  Q.  From  classical  results  regarding  the  representation  of  isotropic 
tensor  functions  (e.g.  [24]),  it  follows  that  N  (a)  can  be  represented  as 

N(a)  =  NiI+N2a+N3a2  (19) 


where  N\,  N2,  and  N$  are  scalar- valued  functions  of  all  stress  invariants.  Third  invariant  dependence 
cannot  be  inferred  without  data  for  which  all  three  principal  stresses  are  distinct,  and  such  data 
are  not  discussed  in  this  paper.  Thus,  we  consider  N(g)  to  be  of  the  form 

N(a)  =  Ni(p,q)I  +  N2(p,q)-  (20) 

<1 


Hence, 


&  =  kT{l 

It  follows  that  for  CTC  loading: 


W(t)  \ 

H(a)J 


Ni(p,q)l+N2(p,q)- 


kTN\(p,q)  = 


kTN2(p,q )  = 


1  _wm 
H(a)  1 

|gj~4i 

l  W(t)\ 

H(a)l 


(21) 


(22) 


(23) 


where  k\  and  £3  denote  the  axial  and  radial  irreversible  strain  rates,  respectively.  Equation  (22) 
indicates  that  N\(p,q )  models  the  irreversible  volumetric  response,  whereas  from  Equation  (23) 
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it  follows  that  N2(p,q)  describes  the  irreversible  shear  response.  We  assume  that  kjN\  can  be 
expressed  as 


kTNi(p,q)  =  <j)(p)  +  il/(p,q) 


(24) 


such  that  0)  =  0.  Thus,  <fi(p)  can  be  determined  using  data  from  the  hydrostatic  compression 
phase  of  CTC  experiments.  Only  the  data  from  the  375  MPa  test  and  the  loading  part  of  the 
450  MPa  cycle  of  the  500  MPa  hydrostatic  experiment  were  utilized.  A  monotonically  increasing 
function  of  pressure  which  tends  asymptotically  to  a  constant  limiting  value  (corresponding  to 
fully  compacted  state) 

<t>(p)  =  <t>  0  + - 

1+exp 

where  </>0  =  — 7.16 x  10-5,  (j>i  =  0.008,  <fi2  =  295. 866 MPa,  and  </>3  =62. 64 MPa,  approximates  the 
data  well.  To  determine  i jj(p,q),  Equation  (22)  along  with  the  data  from  the  deviatoric  phase  of 
CTC  tests  was  used.  Figure  16  shows  the  experimental  variation  of  \ j/  with  the  principal  stress 
difference  (i.e.  q)  obtained  from  various  CTC  tests.  Positive  \ j/  values  correspond  to  compressible 
response,  whereas  negative  ijj  values  indicate  dilatant  behavior.  Note  that  for  unconfined  loading 
the  response  is  largely  dilatant  with  only  very  little  compressibility  at  low  stress  levels.  Based  on 
these  data,  il/(p,q)  was  approximated  with  the  function: 


P~<t>  2 


<A( p,q)  = 


i/'Ci O’-?)  for  q^y(p,q) 

•A c2(P >?)  for  y(p,q)<q^ot(p,q) 
•Ad (/>><?)  for  cc(p,q)<q 


(26) 
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Figure  16.  Experimental  variation  of  kj N\  deviatoric  as  a  function  of  the  principal  stress  difference  q  as 

obtained  from  different  CTC  tests. 
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In  the  above  expression,  y(p,  q)  =  y0ay3l  +y2  with  y0  =  87,  yx  =  0.7,  and  y2  =  1.8,  whereas  oc(p,  q) 
describes  the  compressibility /dilatancy  boundary.  A  decaying  exponential  well  approximates  this 
boundary  (see  Figure  17).  Thus, 


cc(p,q)  =  (Xo  +  (Xi 


^1  —  exp 


(27) 


with  ao  =  3.6MPa,  ot\  =590.05 MPa,  a2  =  0.007,  and  pa  =  IMPa  is  a  dimensionalization  constant. 


The  expressions  of  all  the  other  functions  involved  in  Equation  (26)  are 


<Ac!  (p>  q)  =  P(p,  q)+'l'i  (p>  q) 


1  —  exp 


1  ( q-y(p,<i)\2 

2  V  1A2 Q>>?)  / 


(28) 


exp 


1  /  q-y(p,<i)\2 

2\  h  ) 


1  —  exp 


1  /  ct(p,q)-y(p,q) 

2  \  ^7 


1  /  ct(p,q)-y(p,q) 

2  V  •A  7 


(29) 


and 


•Ad  (/>>«)  = 


OAs^  +  iAs  b)(q~oi(p,q)) 


\-q 


lA9a+lA9f,cr3 
I  +  il/9c  &3 


(30) 


Figure  17.  Compressibility/dilatancy  boundary  for  concrete. 
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with 


P(p,q) 


Po+jhm  I  R 

1  +/^2cr3/ Pa 


•At  (p,q) 


Xhg±h^}/Pa 

i+'At^ 


+{l/ic 


(31) 


/  c3  y  2b 

\l/2(p,q)  =  rl/2a^—J  +^2C 

and  /?0  =  2  x  1(T6,  ft  =  3.3  x  1(T6,  ft =0-4,  ft  =  1.7  x  KT6,  i Ala  =3.3  x  KT6,  iph  =0.004,  i ftc  = 
4.38  xl(T6,  iA2fl  =  1-95,  =0.8,  and  ip3c  =  13,  iA7  =  5GPa,  iA8a  =  -1.2x  KT7,  ifti  =  -3.04x 

10-7,  iA9fl  =60,  \[/c)h  =4.37,  and  i p9c  =0.0043.  The  only  constitutive  function  yet  to  be  determined 
is  krN2(p,q ),  which  governs  the  deviatoric  irreversible  stress-strain  response  (see  Equations 
(23)).  Based  on  unconfined  data  and  CTC  test  data  at  03  =  50,  150,  and  300 MPa,  respectively, 
kTN2(p,q)  was  approximated  by 


kN2(p,q)  =  £  0  — 

,Pa 


1.5 


+  1  —  exp 


q/Pa-<l 

*=3 


u 


(--Qi 

\Pa 


q/Pa-£i 

1  —  (i7/  Pa)/ ^2 


(32) 


All  the  coefficients  involved  in  Equation  (32)  are  functions  of  the  confining  pressure,  i.e. 


03 

Pa 


-  £0  8  +  £oh  ln 


2  (2) 

,(=o 


(33) 


^i(cr3)  =  exp 


—  exp 


°3/Pa-£i  ln(ln(2))-^ 


tic 


and 


^2(^3) : 


+£2h°3/ Pa 


1  +  ^3 /Pa 

In  Equation  (32),  U  is  the  step  function  defined  as 

fl  forv^O 


U(x)  = 


while 


and 


41)(<73)-^0a  + 


42)  (P> «)  =  ?o(,  +  ^0/(1  -exp 


0  forv<0 


^0* 


(l  +  ((<73-^oc)/^)2) 

^3-^oc-^ln(l-V2/2)-^ 


-|\  2 


(34) 

(35) 

(36) 

(37) 

(38) 


The  coefficients  in  the  above  equations  are  £ofl  =  —  2.14x  10-5  ,  ^=2.5xl0-5,  Co,  =  200, 
and  |0d  =  641.51,  £0e  =3.7 x  1CT9,  f0/  =  -2.14x  KT5,  £0g=612.77,  and  <^=441.12,  whereas 
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Zia  =  414.1,  ^  =67.83,  and  ^=61.03,^=61.624,  &*=2.9,  £2c  =0.00177,  and  &  = 
1881464.  A  comparison  between  calculated  and  experimental  kjN2(p,q)  corresponding  to 
unconfined  loading  and  (73  =  150  and  375  MPa,  respectively,  is  shown  in  Figure  18. 

Figure  19  shows  a  comparison  between  calculated  and  experimental  results  for  the  uniaxial 
compression  case.  Note  that  the  model  well  describes  the  overall  behavior.  The  volumetric  response 
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Figure  18.  Comparison  between  experimental  and  theoretical  variations  of  kN2  as  a  function  of  the 
principal  stress  difference  for  <73=0,  150,  and  375  MPa. 
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Figure  19.  Comparison  between  simulated  and  experimental  stress-strain 
responses  under  uniaxial  compression. 
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as  well  as  the  transition  from  compressibility  to  dilatancy  is  particularly  well  captured.  Simulation 
of  the  material  response  under  hydrostatic  compression  is  presented  in  Figure  20.  A  good  correlation 
between  experimental  and  simulation  results  is  obtained.  Finally,  we  present  a  comparison  between 
simulation  and  experimental  results  for  triaxial  compression  at  a  confining  pressure  of  375  MPa. 
Note  that  both  regimes  of  the  volumetric  behavior  (i.e.  compressibility  and  dilatancy)  are  accurately 
simulated  (see  Figure  21). 


Figure  20.  Comparison  between  theoretical  and  experimental  responses  under  hydrostatic  compression. 


Figure  21.  Comparison  between  theoretical  and  experimental  responses  under  triaxial  compression  at  a 

confining  pressure  of  375  MPa. 
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4.  CONCLUSIONS  AND  FINAL  REMARKS 

An  experimental  investigation  into  the  effects  of  high  confinement  on  the  deformation  and  strength 
of  a  concrete  material  was  conducted.  By  performing  CTC  tests  with  several  load-creep-unload  and 
reload  cycles,  we  have  determined  the  elastic  parameters  with  high  accuracy.  The  elastic  parameters 
were  found  to  depend  on  the  stress  state.  Comparison  between  the  mean  stress  vs  volumetric  strain 
data  for  the  deviatoric  phase  of  triaxial  compression  tests  and  the  data  obtained  under  hydrostatic 
compression  led  to  the  following  remarks:  (a)  for  the  same  level  of  the  mean  stress,  the  compaction 
(void  reduction)  under  confinement  is  much  higher  than  under  hydrostatic  compression  and  (b) 
the  maximum  compaction  that  can  be  achieved  depends  on  the  confining  pressure.  It  appears  that 
shear  influences  the  compaction  mechanisms.  Similar  conclusions  have  been  reported  by  Cazacu 
et  al.  [19]  for  alumina  powder  and  recently  by  Burlion  et  al  [25]  for  a  very  fine  aggregate  concrete. 
As  a  consequence,  the  deviatoric  and  hydrostatic  responses  cannot  be  considered  to  be  decoupled. 
In  this  paper,  we  focused  on  modeling  the  effect  of  time  (short-term  creep)  on  the  plastic  response 
at  very  high  confinement  pressures.  An  elastic/ viscoplastic  modeling  approach  was  adopted.  The 
expressions  of  the  creep- stabilization  boundary,  strain-rate  orientation  tensor  (that  governs  the 
evolution  of  the  irreversible  strain)  were  determined  directly  from  the  experimental  data.  We 
compared  the  model  predictions  for  standard  compression  tests  with  data.  Good  agreement  between 
simulated  and  experimental  response  over  the  full  range  of  confining  pressures  (up  to  0.5  GPa) 
was  obtained.  Data  for  non-monotonic  loading  in  this  high-pressure  regime  were  not  available. 
The  validity  of  the  model  for  non-monotonic  loadings  needs  to  be  further  investigated. 


APPENDIX:  SUMMARY  OF  THE  MODEL  PARAMETERS  AND 
CALIBRATION  PROCEDURE 


Elastic  behavior 
Young's  modulus  E 

•  Determined  from  the  slopes  of  the  quasi-linear  parts  of  the  unloading-reloading  cycles 
performed  at  the  end  of  each  of  the  creep  steps  in  the  CTC  tests  at  03  =  300  and  450 MPa, 
respectively. 

•  Parameters  involved  in  the  law  of  evolution  of  E  with  mean  stress  p  (see  Equations  (6)): 
E°°,  b ,  and  c. 

The  Poisson  coefficient  v 

Hydrostatic  yield  function  Hh(p): 

•  Calibrated  by  approximating  the  variation  with  mean  stress  of  the  irreversible  stress  work  at 
creep  stabilization  corresponding  to  the  hydrostatic  phase  of  CTC  tests  at  150,  300,  450,  and 
500  MPa  confinement,  respectively. 

•  Parameters:  an,  bh ,  Ch,  dh ,  and  e\2  (see  Equations  (11)). 

Deviatoric  yield  function  Hd(p,q)\ 

•  Calibrated  by  approximating  the  variation  with  the  principal  stress  difference  q ,  of  Wv> 
(irreversible  work  associated  with  irreversible  shape  change)  corresponding  to  the  deviatoric 
phase  of  CTC  tests  at  150,  300,  450,  and  500  MPa  confinement,  respectively. 

•  Parameters:  doa ,  dob ,  d\a,  d\h ,  d2a ,  d2b,  d2c ,  dsa,  and  d^b  (see  Equations  (13)— (15)). 
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Irreversible  volumetric  response :  N\(p,q)  (see  Equations  (24)) 

•  Its  hydrostatic  component,  denoted  as  <fi(p),  is  calibrated  by  approximating  the  variation  with 
mean  stress  of  the  experimental  irreversible  volumetric  stain  corresponding  to  the  hydro¬ 
static  compression  phase  of  CTC  experiments  at  450  MPa  cycle  and  500  MPa  confinement, 
respectively.  Parameters:  </>0,  </q,  02,  and  03  (see  (25)). 

•  Its  deviatoric  component,  denoted  as  0(p,g),  is  calibrated  by  approximating  the  variation 
with  q  of  the  experimental  irreversible  volumetric  stain  corresponding  to  the  deviatoric  phase 
of  CTC  experiments  at  150,  300,  450  MPa  cycle  and  500  MPa  confinement,  respectively. 

•  Parameters:  y0>  7u  72’  ao>  oci,  and  0C2  (see  (26)-(27)). 

Irreversible  shear  response :  N2(p,q)  (see  Equations  (23)) 

•  Calibrated  by  approximating  the  variation  with  q  of  the  experimental  irreversible  deviatoric 
stain  corresponding  to  the  deviatoric  phase  of  CTC  experiments  at  150,  300,  450  and  500  MPa 
confinement. 

•  Parameters:  Q)a ,  Q,h,  Zoc  and  £0e,  £0/,  £0g,  Co,, ,  £ia,  £ib,  £ic,  and  £3. 
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